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ABSTRACT 



Aims. We investigate the molecular gas properties of a sample of 23 galaxies in order to find and test chemical signatures of galaxy 
, evolution and to compare them to IR evolutionary tracers. 

Methods. Observation at 3 mm wavelengths were obtained with the EMIR broadband receiver, mounted on the IRAM 30 m telescope 
on Pico Veleta, Spain. We compare the emission of the main molecular species with existing models of chemical evolution by means 
of line intensity ratios diagrams and principal component analysis. 
f-\ , Results. We detect molecular emission in 19 galaxies in two 8 GHz- wide bands centred at 88 and 112 GHz. The main detected 

^ , transitions are the 7=1-0 lines of CO, '''CO, HCN, HNC, HCO+, CN, and CoH. We also detect HC,N 7=10-9 in the galaxies 

^ ■ IRAS 17208, IC 860, NGC 4418, NGC 7771, and NGC 1068. The only HCjN detections are in objects with HC0+/HCN<1 and warm 

' IRAS colours. Galaxies with the highest HCjN/HCN ratios have warm IRAS colours (60/100 ytim>0.8). The brightest HCjN emission 

is found in IC 860, where we also detect the molecule in its vibrationally excited state. We find low HNC/HCN line ratios (<0.5), that 
cannot be explained by existing PDR or XDR chemical models. The intensities of HCO+ and HNC appear anti-correlated, because 
galaxies with low HCO+/HCN intensity ratios have high HNC/HCN. No correlation is found between the HNC/HCN line ratio and 
^ ' dust temperature. All HNC-bright objects are either luminous IR galaxies (LIRG) or Seyferts. Galaxies with bright polycyclic aromatic 

, hydrocarbons (PAH) emission show low HNC/HCO^ ratios. The CO/'^^CO ratio is positively correlated with the dust temperature and 

OA ■ is generally higher than in our galaxy. The emission of CN and C"*0 is correlated. 

Conclusions. Bright HC3N emission in HCO^-faint objects may imply that these are not dominated by X-ray chemistry. Thus the 
OA ' HCN/HCO^ line ratio is not, by itself, a reliable tracer of XDRs. Bright HC3N and faint HCO^ could be signatures of embedded star- 

formation, instead of AGN activity. Mechanical heating caused by supernova explosions may be responsible for the low HNC/HCN 
and high HCO^/HCN ratios in some starbursts. We cannot exclude, however, that the discussed trends are largely caused by optical 
depth effects or excitation. Chemical models alone cannot explain all properties of the observed molecular emission. Better constraints 
to the gas spacial distribution and excitation are needed to distinguish abundance and excitation effects. 

Key words, galaxies: evolution — galaxies: starburst — galaxies: active — radio lines: ISM — ISM: molecules 
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X. 

^1. Introduction AGN and starburst are still not well understood and must be 
■ - - ' further explored. At high redshift, LlRGs dominate the cos- 
Luminous infrared galaxies (LlRGs) radiate most of their lu- j^ic infrared background and, by assuming that they ai-e pow- 
mmosity (Lir >10 Lq) as dust thermal emission in the gj-ed by starburts, we can use these galaxies to trace the dust- 
infrai-ed and have been studied at almost all wavelenghts obscured star-formation rate, th e dust c ontent and the metal- 
d Sanders&Mirabell 1 19961) . However, the nature of the power jj^ity early Universe (Bar ger et al.l [l999). The most well 
source is still unclear when the inner region of the LIRG is known techniques to distinguish between AGN and star-powered 
obscured by dust. The high central IR surface brightness im- galax ies rely on the observation of emission lines in the optical 
plies that this power source can be either an embedded com- (g IVeilleux & Osterbrockll987h . Over the last decade several 
pact starburst or an enshrouded AGN - or a combination of Q^her diagnostic diagrams, based on IR spectra, have been pro- 
both. The evolution of the activity and the connection between ^^^^^ quantify the contribution of star-formation and AGN 

activit y to the infrared lum inosities of LlRGs ( Genz el et al.l 

* F.C. wishes to thank the EU ESTRELA programme for support. 1998t iLutz et anil998t ISpoon et al. 2007) . In the mm and sub- 

** S.A. wishes to thank the Swedish Research Council for grant sup- mm there have been attempts as well to classify the activ- 
port. 
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ity of galaxies via diagnostic diagrams, of whi ch the most 
well k nown is the HCN/HCO+ line ratio plot of iKohno et all 
(1200 lb and Jmanishi et al. (2004). A multi-transition study of 
the HCN/HCO^ ratio in Seyfert and starburst galaxies has also 
been reported by IKrips et al.l (l2008l) . These authors find an un- 
derluminosity of HCO^ in some AGN-dominated cores, which 
they suggested to be owing to an underabundance of HCO^ 
caused by the X-ray-d ominated chemistry induced by the AGN 
dMalonev et alj IT996I) . This latter interpretation has been dis- 
puted since it does not agree with recent models of X-ray- 
dominated regions (XDRs) (Meiierink& Spaans 2005), which 
instead show an enhancement of HCO^ abundances caused by 
an increase in ionization. Other interpretations in terms of star- 
burst evolution h a ve been put forward as an alternative (e.g., 
iBaan et all |2008|) . iGracia-Carpio et alj (l2006h also found that 
ULIRGs (LiR >10'^ Lo) in general seem to have lower HCO^ 
1-0 luminosities with regards to HCN - compared with more 
moderate LIRGs. The HCN/CN 1-0 and HCN/HNC 1-0 line ra- 
tios have been used as w ell to help interpret a galaxy's pos ition in 
an evolutionary scheme dAalto et al.'2002':l Baan et alj200 8). For 
the brightest nearby gala xies, a sch eme based on the HNCO/CS 
ratio has been proposed dMartm et al.ll200 9). Moreover, promi- 
nent sources such as NGC 253 and IC 342 allow the detec- 
tion of rarer species that can be used in turn to help identifying 
and even resolving the dominant activity (e.g., iMeier & Turnen 
[2005; Martin et al. 2006). 

Line ratios need to be very accurately measured if we want them 
to be sensitive tracers of molecular properties. The new EMIR 
receiver, mounted in 2009 on the IRAM 30 m telescope in Spain, 
offers the opportunity of achieving this high accuracy. The avail- 
able bandwidth of nearly 8 GHz at 3 mm allows us to fit many 
key molecular lines in the same band, therefore eliminating the 
uncertainties owing to relative calibration and pointing errors 
that affect most single-dish observations. We present a survey 
of molecular lines for a sample of 23 galaxies observed with 
EMIR in the period June-November 2009. In section |3] we re- 
port the details about the observations and source selection. In 
Sect. ID we present the results. In Sect. |5]the line ratio diagrams 
for HCN, HNC, HCO+, and HC3N are presented and discussed. 
The molecular emission is also compared with far-infrared (FIR) 
colours and polycyclic aromatic hydrocarbons (PAH) emission. 
In Sect. |6]we present our conclusions and an outlook. The spec- 
tra of all observed sources and the tables summarizing the line 
parameters of the detected species can be found in the Appendix. 

2. Source selection 

The first diagnostic diagram to take into account the effects of 
str ong obscuration of the nuclear power source was presented 
bv lSpoon etan(l2007h . using the equivalent width of the 6.2 yum 
PAH emission feature and the strength of the 9.7 jjm silicate ab- 
sorption (see Fig.[TJ. 

Based on the position in the diagram, galaxies are put into nine 
classes, ranging from continuum-dominated AGN hot dust spec- 
tra (lA) to PAH-dominated starburst spectra (IC) to absorption- 
dominated spectra of deeply obscured galactic nuclei (3A). 
[Spoon et al. (2007) find that galaxies are systematically dis- 
tributed along two distinct branches: one horizontal sequence at 
low silicate depths, ranging from AGN to starburst-dominated 
spectra, and one diagonal sequence at higher silicate strength, 
ranging from obscured nuclei to pure-starburst objects. The sep- 
aration into two branches likely reflects fundamental differences 
in the dust geometry in the two sets of sources. Spectra of lu- 
minous infrared galaxies are found along the full length of both 
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PAH equivalent width (6.2 |im) [|im] 



Fig. 1. Position of the o bserved galaxies on the mid-infrared 
diagnostic diagram from ISpoon et aP (l2007l) . Galaxies are di- 
vided into nine classes (1A-3C), depending on their spectral 
properties and position on the graph. These classes range from 
PAH-dominated spectra (IC) to continuum-dominated spec- 
tra with faint PAH emission (lA) and obscured galaxies with 
deep silicate absorption (3A). Galaxy types are distinguished by 
their plotting symbol: Squares: LIRGs. Stars: Starburst galax- 
ies. Triangles: Seyfert galaxies. The numbers indicate differ- 
ent galaxies, as explained in Table |2] The shadow ed area cor- 
responds to the two main sequences discussed by I Spoon et"an 
(2007). An increase of the PAH equivalent width indicate an in- 
crease of the starburst contribution to the emission, while AGN 
spectra are continuum-dominated. See Section|2]for further dis- 
cussion. 

branches, reflecting the diverse nature of the LIRG family. In 
this work, the term LIRG broadly refers to objects which emit 
most of their energy in the IR. Thus a galaxy classified as LIRG 
may be a starburst, an AGN, or both, but the dust obscuration 
hinders a clear classification. We do not explicitly distinguish 
between LIRGs (Ljr >10" Lq) and ULIRGs (Ljr >10^^ Lq). 
However, IR luminosities are reported in Table |2] for reference. 
Objects with both LIRG and AGN signatures are labelled as A, L 
in Tablel and appear as AGNs (triangles) in the graphs. 
ISpoon et al. (.2007.) interpret the observed distribution as a pos- 
sible evolutionary effect, with sources moving from the diagonal 
to the horizontal branch as the dust distribution evolves from 
a uniform to a clumpy geometry. The underpopulated 2A class 
implies either that this transition for LIRGs is very rapid, or that 
LIRGs mostly evolve into unobscured starbursts. 
We aim to compare this mid-IR evolution scheme with mm 
molecular observations. Our targets were mainly selected from 
the sample of Spoon et al. (2007), plus a few interesting objects 
for which we had PAH and silicate mid-IR data. The source se- 
lection criteria were the following: 

- Uniform coverag e of the most signific ant classes in the diag- 
nostic diagram o ASpoon et al.\ i200'A) . 

- Uniform representation of different galaxy types ( Seyferts, 
starbursts, LIRGs), 

- Source visibility at IRAM site. 

The resulting sample is composed of 23 galaxies, whose 
main properties are listed in Tab |2] All targets were observed 
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„ I „ , „ HCO+(I-0) HNC(l-O) HC,N(10-9) C2H(l-0) CO(l-O) t /t /t \ t m 

Galaxy RA Dec VHei.o hcn(I-O) hcn(I-O) hcn(I-O) hcn(I-O) ocod-O) LogCLlR/Lo) Class Type Num. 

[h:m: s] [°:':"] [km/s] 



IRAS 17208 


17:23:21.9 
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12834 


IC 860 


13:15:03.5 


+24:37:08 


3347 


Mrk231 


12:56:14.2 


+56:52:25 
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NGC 1614 


04:33:59.8 


-08:34:44 


4778 


NGC 3079 


10:01:57.8 


+55:40:47 


1116 


NGC 4194 


12:14:09.5 


+54:31:37 


2501 


NGC 4388 


12:25:46.7 


+ 12:39:44 


2524 


NGC 4418 


12:26:54.6 


-00:52:39 


2110 


NGC 6090 


16:11:40.7 


+52:27:24 


8785 


NGC 6240 


16:52:58.9 


+02:24:03 


7339 


NGC 7469 


23:03:15.6 


+08:52:26 


4892 


NGC 7771 


23:51:24.9 


+20:06:43 


4277 


NGC 660 


01:43:02.4 


+ 13:38:42 


850 


NGC 3556 


11:11:31.0 


+55:40:27 


699 


NGC 1068 


02:42:40.7 


-00:00:48 


1137 


NGC 7674 


23:27:56.7 


+08:46:45 


8671 


UGC 2866 


03:50:14.9 


+70:05:41 


1232 


UGC5101 


09:35:51.6 


+61:21:11 


11802 


NGC 2273 


06:50:08.6 


+60:50:45 


1840 


Arp220 


15:34:57.2 


+23:30:09 


5382 


IRAS 15250 


15:26:59.4 


+35:58:38 


16535 


NGC 1 140 


02:54:33.6 


-10:01:40 


1501 


NGC 1056 


02:42:48.3 


+28:34:27 


1545 


NGC 1377 


03:36:39.1 


-20:54:07 


1792 
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0.77(0.14) 


0.31(0.07) 


0.42(0.10) 


0.62(0.16) 


0.62(0.16) 


0.42(0.13) 


1.04(0.29) 


0.56(0.08) 


0.38(0.06) 


<0.07 


0.29(0.12) 


1.83(0.37) 


0.33(0.15) 


<0.35 


1.00(0.33) 


1.12(0.11) 


0.27(0.05) 


<0.06 


0.54(0.09) 


1.32(0.35) 


0.53(0.23) 


<0.20 


0.98(0.36) 


1.38(0.40) 


0.62(0.24) 


<0.35 


0.81(0.52) 


0.59(0.10) 


0.47(0.09) 


0.37(0.09) 


0.64(0.18) 


1.67(0.44) 


0.25(0.17) 


<0.18 


0.46(0.20) 


1.63(0.14) 


0.20(0.06) 


<0.09 


0.36(0.08) 


1.12(0.11) 


0.55(0.08) 


<0.07 


0.75(0.24) 


0.94(0.09) 


0.44(0.07) 


0.06(0.04) 


0.32(0.06) 


1.04(0.09) 
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<0.27 


1.56(0.37) 


0.67(0.03) 


0.40(0.03) 


0.04(0.02) 


0.37(0.09) 


1.46(0.20) 


0.51(0.13) 


<0.14 


0.88(0.17) 


0.36(0.20) 


0.82(0.28) 


<0.28 


0.89(0.37) 


1.05(0.37) 


1.09(0.38) 


<0.89 


<0.89 


0.47(0.07) 


0.49(0.12) 


0.19(0.06) 







12.35 


IB 


L 


1 


18.96(3.83) 


11.14 


2B 


L 


2 




12.37 


lA 


A,L 


3 
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11.43 


IC 


S 


4 
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10.65 


2C 


A 


5 
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IC 


S 
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9.66 


2A 


A 


7 




11.00 


3A 


L 
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IC 


L 


9 


28.79(2.99) 
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2B 


L 


10 
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11.41 


IB 


A 
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13.61(0.52) 
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L 


12 


16.45(0.49) 
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2C 


S 


13 


12.50(0.31) 


10.00 




s 
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10.89 


lA 


A 
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14.46(1.18) 
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A 
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S 
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2B 


A,L 
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IB 


A 


19 
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3B 


L 


20 
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3A 


L 
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IC 


S 


22 
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S 
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3A 


o 
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Table 2. Line ratios and general properties of the observed galaxies. The mid-IR Class column refers to the classification by 
[Spoon et al. (2007), based on the PAH equivalent width and silicate absorption in the mid-infrared. The column Type summarizes 
the properties of the object, distinguishing between starbursts (5'), Seyferts (A), and LIRGs (L). The galaxy NGC 1377 is classified 
as obscured (O), because it emits mostly in the IR, but its luminosity is not high enough to be classified as LIRG. For a discussion 
about source types see Section|2] In Figs. [Till] and |5] the sources are labelled with the numbers reported in column A^mot. 



in the 88 GHz band, but because of time constraints, only 12 
were observed in the 112 GHz band. 



3. Observations 

3.1. Line ratios at 88 GHz 

The observations were obtained in June-November 2009 with 
the IRAM 30 m telescope on Pico Veleta, Spain. The 8 GHz 
band of the EMIR receiver was centred on two different tunings, 
at 88.675 and 112.15 GHz. These frequencies were chosen in or- 
der to accommodate as many potentially strong lines as possible 
in the same band. A list of the transitions of the most important 
molecular gas tracers is given in Table [1] together with beam ef- 
ficiencies (77mb) and beam sizes (HPBW) for the two observed 
bands. Observations were performed in dual polarization, each 
covering 8 GHz, which nicely fitted the 4x4 GHz backend bot- 
tleneck. The E090 frontend was connected to the low-resolution 
WILMA autocorrelator, which is capable of processing the large 
input bandwidth with a resolution of about 2 MHz (^7 and 5 
km s"' at 88 GHz and 112 GHz, respectively). Observations 
were performed in wobbler switching mode, with a throw of 60- 
120" (depending on source size), in order to maximize baseline 
quality. 

The pointing model was checked against bright, nearby calibra- 
tors for every source, and every two hours for long integrations. 



Calibration scans on the standard two load system were taken 
every 5 minutes. The observing conditions were optimal, with 
characteristic system temperatures of 100 and 150 K for the 88 
GHz and 112 GHz observations, respectively. This resulted in a 
rms noise per channel of roughly 0.6 mK at 60 km s ' resolu- 
tion, for 2 hours of on-source observing time. 
Data were reduced with the CLASqJ software. A first order 
baseline was removed from all spectra, which are shown in 
Appendix ini The intensity scale in all spectra is in T^, which 
is related to the source brightness temperature as 



t: el + el 

rimh et 



(1) 



with 6*5 and e^, the angular sizes of source and beam, respec- 
tively. We assume a constant main beam efficiency across both 
observed bands. The error on the main beam brightness temper- 
ature estimate introduced by this assumption is of the order of 
one percent. 

4. Results 

Rest frequencies were taken from the NIST database 
Recommended Rest Frequencies for Observed Interstellar 



http://iram.fr/IRAMFR/GILDAS/ 
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Fig. 2. Spectra observed with EMIR in IC 860 at 88 and 1 12 GHz. The intensity scale is in T* , not corrected for main beam efficiency. 
The region marked with (1/10), around the CO 1-0 line, has been scaled down by a factor 10. The main molecular transitions are 
labelled regardless of line detection. The C2H 3/2 and 1/2 labels mark the limits of the C2H multiplet at 87 GHz. Transitions of 
vibrationally excited HC3N are labelled as le and If. The frequency scale is the observed frequency, not corrected for redshift. 



Molecular Microwave Transition^ Our line identification takes 
into account the distortion of the velocity scale caused by the 
large observed band width. For a discussion of this effect, see 
lGordonetalJ(ll992h . 

An example of an EMIR spectrum is shown in Fig. |2l while all 
the observed spectra are shown in Appendix iDl Line intensities 
were extracted by means of Gaussian fitting and are reported 
along with other line parameters in Appendix|El 
Integrated intensities of the most relevant species were combined 
to form line intensity ratios, which are listed in Table |2] These 
were used to produce the diagrams shown in Figs. |4] and |5] We 
choose to compare intensities of transitions inside the same fre- 
quency band to maximize the accuracy of the derived line ratios. 
The diagrams derived from the two observed bands are described 
in sections [3n and l4TT] A principal component analysis of the 88 
GHz dataset is discussed in section l4!2l 

In plots a, b, and c of Fig. [H we report the line ratios between 
the first transitions (7=1-0) of HCN, HNC and HCO+ for all de- 
tected galaxies. The different symbols refer to different galaxy 
types. We note that all non-compact starburst galaxies (stars) 
are HCO+ -luminous, with line ratios HCO+/HCN > 1 (plots a, 
b). Luminous infrared galaxies (squares) and Seyferts (triangles) 
do not show any strong trend. However, the majority of active 



galactic nuclei do reside at low HCO^/HCN values compared 
with starbursts. 

The HC3N/HCN line ratio is also reported in the graphs. Circles 
are drawn around sources where HC3N 7=10-9 has been de- 
tected, the diameter of the circle being proportional to the 
HC3N/HCN ratio. Evidently that all the HC3N detections have 
HCO^/HCN< 1 and HNC/HCN>0.4. 

Plots a and c show strong correlations. As HNC/HCO^ in- 
creases, galaxies move on the graphs towards higher HNC/HCN 
and HCO^/HCN ratios. An inverse correlation between HNC 
and HCO^ line intensities also emerges from plot b. Galaxies 
at HCO^/HCN < 1 present on average a 30 % increase in 
HNC/HCN, compared with those at higher HCO+/HCN ratios. 
If we exclude the Seyfert NGC 2273, the higher end of the HNC 
luminosity distribution (see graphs b and c in Fig. lU is com- 
posed mainly of luminous infrared galaxies. 
The general picture emerging from the graphs is that starbursts 
are characterized by faint HNC and bright HCO^ emission 
(compared with HCN), while LIRGs mostly occupy the oppo- 
site end of the HNC-HCO^ coiTelation, with high HNC and low 
HCO^ intensities. All HNC-bright galaxies are either LIRGs or 
Seyferts. Most of the HC3N detections are LIRGs, the only ex- 
ception being the Seyfert NGC 1068. 



^ http://physics.nist.gov/PhysRefData/Micro/Html/contents.html 
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Table 1. Properties of the two observed bands. Rest frequencies 
of the brightest Hnes in each band are shown with the beam sizes 
(HPBW) and main beam efficiencies (r\^h)- Energies of upper 
level (Eup) and critical densities at 20 K (n^ were taken from the 
Leiden atomic and molecular database. 



HNC/HCN HCO+/HCN C2H/HCN Silicate PAH 25/100 pm 



HNC/HCN 


1.00 












HCO+/HCN 


-0.43 


1.00 










C2H/HCN 


0.46 


0.09 


1.00 








Silicate 


-0.11 


-0.37 


0.03 


1.00 






PAH 


-0.41 
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0.19 
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1.00 




25/100 
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-0.10 
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0.03 
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1.00 



Table 3. Correlation matrix 



4 J. Une ratios at 1 12 GHz 

The main detections in the 112 GHz band are the J-l-0 emission 
hnes of CO, '^CO and C'^O, and the spin doublet of CN 1-0 
J -3/2-1/2 and 7=1/2-1/2. These lines are detected in most of 
the galaxies in our sample, and their integrated intensities are 
compared in the plots in Fig.H] 

A particularly rich chemistry is found in IC 860, whose spectrum 
shows bright emission lines of methanol and HC3N 7=12-11. 
This detection of HC3N is the only one in the 1 12 GHz sample 
and will be further discussed in section l5l8] 

4.2. Principal component analysis 

Line-ratio diagrams as the ones shown in sections 13.11 and 14.11 
represent the standard framework for interpreting molecular 
data and can be easily compared with previous studi es (e.g., 
iBaan et al.l l2008: 'Loen en e"tani2008l: ISaan et alJ l20ltf). as dis- 
cussed in section l43] In this approach, the properties of different 
galaxy types are compared two at a time to find the observables 
that best characterize different physical environments. In gen- 
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PC 2 


PC 3 


PC 4 


PC 5 


PC 6 


Variance /c 
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16 




3 




HNC/HCN 


0.03 


-0.19 


0.59 


-0.32 


0.02 


-0.71 


HCO+/HCN 


-0.24 


0.17 


-0.26 


-0.18 


-0.9 


-0.21 


C2H/HCN 


-0.02 


0.10 


0.26 


-0.79 


-0.01 


0.54 


Silicate 


0.83 


0.52 


-0.06 


-0.06 


-0.07 


-0.12 


PAH 


-0.46 


0.68 


-0.20 


-0.17 


0.42 


-0.28 


25/100^;m 


0.17 


-0.43 


-0.69 


-0.45 


0.22 


-0.24 



Table 4. Projection of the principal components onto the base 
of observables. For each PC, its contribution to the total variance 
in the data is also shown. 



eral, however, the observed quantities are not independent and 
different environments are best described by a combination of 
observables. As the number of observables increases, it becomes 
more and more difficult to interpret multidimensional datasets by 
means of 2D sections (i.e., line-ratio diagrams). 
To address this complexity, we applied a principal component 
analysis (PCA) to our data. This is commonly used to reduce the 
dimensionality of a dataset and is very effective in finding hidden 
trends that could otherwise be buried in the noise. For an appli- 
cation of a PCA to molecul ar map s, see, e.g., Ungerechts et aT| 
(I1997I) and lMeier & TurneJ(l2005l) . 

Each galaxy in our sample is described by a set of n observed 
quantities (e.g. line ratios, IR properties) that represent an ini- 
tial base of vectors. The PCA algorithm first computes the co- 
variance matrix of the data along the n directions and finds its 
eigenvalues and eigenvectors. The eigenvalues are then sorted in 
descending order, and the corresponding eigenvectors labelled 
as principal component (PC) 1, 2, etc... 

The first principal component (PC 1) is thus the linear combina- 
tion of the initial galaxy properties along which the dispersion 
of the data is maximum. PC 2 is the vector, perpendicular to PC 
1, which has the second highest dispersion, and so on for all n 
PCs. The values of the observed properties for each galaxy are 
then projected on the new base of PCs. 

In our analysis we choose as a base of observables the line in- 
tensity ratios HNC/HCN, HCO+/HCN, and C2H/HCN and the 
IR properties silicate absorption, PAH EW and the ratio of the 
IRAS fluxes at 25 and 100 yum. Because our algorithm cannot 
deal with upper limits, we do not include faint lines, as HC3N. 
In order to maximize the number of galaxies, we limit our anal- 
ysis to the 88 GHz band, since only a fraction of our sample was 
observed at 1 12 GHz. 

4.2.1 . Results of the PC analysis 

The results of the PCA algorithm are shown in Table |4] Here the 
projection of the PCs onto the original base are reported along 
with their contribution to the total dispersion in the dataset. Most 
of the information is contained in the first three PCs, which ac- 
count for about 90 % of the total variance. In the following anal- 
ysis we will thus focus on PC 1,2,3. 

In Fig. [3] right we show the projection of the PCs onto the origi- 
nal base of observables, while in Fig. [3] left we plot the position 
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Fig. 3. Results from the principal component analysis. Left: Position of the observed galaxies in the PC base [PC\,PC2,PC3]. 
Different symbols refer to starbursts {stars), Seyferts {triangles) and LIRGs {squares). Galaxies with both AGN and LIRG properties 
are p lotted as triangles. N umbers indicate different galaxies, as explained in Table|2] Values for the PAH equivalent widths are taken 
from lSpoon et alj ( 2007h . Infrared fluxes were obtained from the NASA/IPAC Extragalactic Database. Right: Projection of the PC 
components onto the original base of observables. We consider only the first three PCs, since they account for 90 % of the variance 
and thus contain most of the information in the dataset. 



of the galaxies in the PC base. 

The main contribution to the dispersion of the data comes from 
PAH EW and silicate absorption, which have the highest projec- 
tions on PC 1 . The opposite signs mean that the dispersion along 
PCI is mainly caused by galaxies that have low PAH EW and 
high silicate absorption, or vice versa. The first PC is therefore 
dominated by the diagonal sequence on the mid-IR diagram of 
Fig. [U Galaxies with high values of PC 1 will have high silicate 
absorption and low PAH EW. This is evident when comparing 
Figs.|I]and[3] left, where the galaxies NGC 4418 and NGC 6090 
(8 and 9 on the graphs) lie at the opposite ends of both PC 1 and 
the mid-IR sequence. 

Among the molecular line intensity ratios the largest contri- 
bution to PC 1 is given by HCO+/HCN. From Fig. [3] left, 
HCO+/HCN appears to be positively correlated with PAH EW 
(same sign of the projection along PC 1). This agrees with 
Figs. IIJ/ and IB. II and generally with the trend that will be dis- 
cussed in Sect. 15.101 

The second principal component (PC 2) is also dominated by IR 
observables, but with diff'erent contributions than the ones de- 
rived for PC 1 . Silicates and PAH projections are now positively 
correlated, and the ratio IRAS 25/100 yum gives a significant con- 
tribution to the dispersion. As can be seen in Fig. [3] left, PC 2 ef- 
ficiently separates AGN s from obscured or starburst- dominated 
sources. All AGNs, with the exception of NGC 3079 (number 
5 on the graph), have negative PC 2 values, which correspond 
roughly to classes lA-lB on the mid-IR diagram of Fig. [T] An 
interesting anti-correlation between PAH EW and IRAS 25/100 
jum flux ratio is also evident in Fig. [3] left. This may be caused 
by higher dust temperatures in systems with a dominant AGN 



component. 

Most of the molecular information is contained in PC 3, 
where a substantial fraction of the total dispersion is along 
the HNC/HCN-HCO+/HCN direction, while the contribution 
by PAH and silicates is almost negligible. The IRAS colour 
(25/100 /im ratio) still plays an important role, with a projection 
along PC 3 (--0.7) slightly higher in absolute value than the one 
for the HNC/HCN line intensity ratio (^0.6). The opposite sign 
of these two projections is caused by a slight anti-correlation of 
the HNC/HCN ratio with dust temperature, which also results 
from the correlation matrix in Table |3] This trend, however, is 
not evident in Fig. |4^. The distribution of galaxies along PC 3 
does not show any obvious correlation with galaxy properties 
and its interpretation is not straightforward. We will further dis- 
cuss the results of the PC analysis in Paper II, when we will have 
additional information from VLA radio observations about, e.g., 
gas surface density and star-formation rate. 



4.3. Comparison with previous observations 

A previous study of the molecular emission of dense gas in lu- 
minous infrared galaxies was reported by Baan et al. (2008). The 
authors analyse data for 117 galaxies, with infrared luminosities 
ranging over about three orders of magnitude. Our diagrams in 
Fig. Scan be directly compared with those of Fig. 9 in lBaan et al.l 
(120081) . 

The correlations between HNC/HCO^ and the two ra tios 
HNC/HCN and HCO+/HCN, shown by iBaan et all dloo^ in 
plots a and c, are confirmed by our observations. Our measure- 
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ment have a much smaller scatter, thanks to the excellent relative 
calibration of the line intensities provided by the EMIR receiver. 
However, we cannot exclude that the smaller scatter is partially 
caused by source select i on effe cts. 

Plot b in .Baan eTal] (l2008h . showing HCO+/HCN versus 
HNC/HCN, however, cannot be reproduced by our data. In par- 
ticular, we did not observe any galaxy with log HNC/HCN<-0.4 
an d log HCO+/HCN <0, while about 25 % of the sources plotted 
bv lBaan et al.l (l2008h are in this range. 

As a result, the correlation between HCO+/HCN and HNC/HCN 
line ratios in our graph b goes in the opposite direction, with low 
HNC/HCN corresponding to high HCO+/HCN. It is not clear 
whether this is caused by selection effects tha t may affect our ob - 
servations, or to the large scatter in the data in lBaan et al.l (l2008h . 
In Table lC.ri we compare line intensities from^ Baan et al.l (l2008l) 
with our EMIR values. In most cases the two datasets agree at 
the 20 % level, with some exceptions. The most striking differ- 
ences are seen for HNC/HCN in NGC 7469 and HCO+/HCN 
in NGC 6240, which vary more than 60 %, passing from values 
higher than unity in one dataset to much lower ratios in the other. 
These discrepancies could be caused by calibration or pointing 
inaccuracies, which affect most narrow band single dish observa- 
tions. The determination of line ratios with EMIR is more robust, 
because these inaccuracies do not affect the relative intensity of 
lines detected in the same band. 



5. Discussion 

5.1. Interpreting line-intensity ratios 

Line-intensity ratios are often interpreted as estimates of the 
abundance ratios of two molecular species. This is true only if 
the emission is optically thin and the two molecules have simi- 
lar excitation properties. Moreover, single dish observations only 
reveal global ratios, averaged on the whole galaxy, and the inter- 
pretation of the results has to take into account the lack of spacial 
information. 

In our study we detect transitions with critical densities (wc) that 
vary by almost four orders of magnitude (see Table [1), and it is 
therefore likely that excitation plays an important role in deter- 
mining the observed ratios. 

Transitions with large ric are efficiently excited in high density 
gas, which is mostly concentrated in the inner parts of galaxies, 
while the emission of lines with lower n^, which can be excited in 
more diffuse gas, is generally more extended. Low-Wc lines will 
thus have a higher beam filling factor, resulting in an increased 
T^. In this context, it is advisable to compare line intensities of 
transitions with similar Wc, since the resulting ratios can be more 
directly related to molecular abundances. 

Also, the excitation temperature of the transitions is an impor- 
tant factor when interpreting single dish intensities. If the emis- 
sion is coming from regions with different gas temperatures, this 
translates into a different population of the low-J levels, which 
directly affects the line intensity ratios. 

Another complication arises when we take into account radia- 
tive excitation. Many of the observed molecules, as e.g. HNC, 
HCN, and HC3N, can in fact be affected by radiative pumping 
via IR vibrat ional modes . Radiative pumping of HNC has been 
proposed by lAalto et aP (12007.) as an explanation of the high 
HNC/HCN ratio observed in LIRGs, while the strong connection 
of H C3N excitation with the IR continuum has been discussed in 
ICostagliola & Aalto (2010). 

Extreme caution should be used therefore when interpreting 
global line intensity ratios as abundance indicators. Molecular 



excitation across our galaxy sample will be further discussed in 
a following paper, where we will include 1 mm counterparts of 
the observed lines. 



5.2. Is HCO+/HCN driven by XDR chemistry? 



The HCO^/HCN J-l-0 line ratio has been proposed as a reli- 
able diagnostic tool to disting uish betwe en AGN- and starburst- 
powe red galaxies (e.g., iKohno et al.li200 It Jmanishi et al. 20041 
I2OO7I) . Interferometric observations by iKohno et al.. (,200 Ih 
reveal that Seyfert galaxies have lower HCO^/HCN ratios 
compared with starbursts. Seyfert galaxies with starburst-like 
HCO^/HCN rati os are generall y interpreted as mixed AGN- 
starburst objects dlmanishi et al.l2 007). In our observations, star- 
bursts and LIRGs have, on average, higher HCO^/HCN ra- 
tio s than AGN , whic h is consistent with the trend observed 
by Kohno et alj (l200lh and llmanishi etalJ ('2007). This is com- 
monly attributed to an enhancement of HCN abundance in the 
X-ray-dominat ed region (XDR, |M alonev et al.ll 1996) surround- 
ing an AGN ("Lepp & Dal garnoi il996i) . Chemical models by 
Meijerink & Spaans (2005) show that an HCO+/HCN line ratio 
lower than one can be observed on the surface of low-density 
(n <10^ cm""*) XDRs, for H2 column densities <10^^ cm"^. 
However, at higher depths inside the molecular cloud, HCO^ 
abundance increases and eventually the cumulative HCO ^/HCN 
column density ratio becomes higher than one. Meiierin k et al.l 
hOOlh report that for N(H2)>10^'' cm"^, the intensity of the 
7=1-0 and J -4—3 transitions emerging from an XDR is higher 
for HCO^ than for HCN by a factor of at least three. This is at- 
tributed to the fact that HCN abundance in an XDR is enhanced 
(compared with HCO^ ) only in a narrow range of ioniza tion 
rate/de nsity (Fig. 3 in ILepp & Dalgarnolll996l) . iMeiierink et^H 
(120071) conclude that a low HCO+/HCN hne ratio is a good 
tracer of UV-dominated regions (PDR) for gas densities >10^ 



cm rather than X-ray-irradiated gas. Thus the observed low 
HCO^/HCN line ratio could be either caused by low-density 
(« <10^ cm.-^) XDRsorderise (n > 10^ cm"^) PDRs. 
Observations bv lBlake et aLld 1987 b in the Orion molecular cloud 
reveal extremely low HCO^/HCN abundance ratios (<10"'') in 
hot cores, i.e. warm and dense gas aro und massive young stars. 
This is supported by calculations by iBavet et aP (12008). who 
find abundances of HCO^ as low as 10"'^ (relative to H2) in 
a wide range of hot core models. The faint HCO^ emission may 
thus be caused by deeply embedded star-formation, instead of 
XDR or PDR chemistry. 

This interpretation is based on the assumption that line ratios are 
directly linked to relative abundances. This may not be true for 
HCN and HCO^, which have different excitation properties. The 
critical density of the 7=1-0 transition of HCO^ is one order of 
magnitude lower than for HCN. As a consequence, HCO^ emis- 
sion can originate from lower density gas, which in general has 
a higher filling factor compared with the dense («h >10^ cm"^) 
component. This implies that the high HCO^/HCN line ratio ob- 
served in starburst galaxies may be caused by different filling 
factors fo r the two molecules instead of by their abundance. 
Although Ixi-ips etal.l (l2008l) suggest that the HCN/HCO+ line 
ratio can often serve as a measure of abundance, we advise cau- 
tion in doing so. Recent work by Sakam oto et al.l (12009 ^ and 
Aalto et al. (2009) show that in luminous, dusty nuclei absorp- 
tion and radiative excitation affect the line emission of these 
molecules on global scales. 
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Fig. 4. Diagrams derived from observed line ratios. Symbols and galaxy labelling are the same as in Fig. [3] Galaxies where HC3N 
was detected are marked with a circle. The diameter of the circle is proportional to the HC3N 10-9/HCN 1-0 line ratio. Error bars 
show l-o" uncertainties. 



5.3. The HNC/HCN ratio 

In the Galaxy, the H NC/HCN observed li ne ratio ranges from 
1/100 in hot cores ([Schilke et al.' '1992) to values as high 
as 4 in dark clouds (Hirota et al. 1998). The abundance of 
the HNC molecule decreases with increasing gas temperature. 
[Hirota et al. ( 1998) suggest that this may be owing to the temper- 
ature dependence of neutral-neutral reactions, which, for tem- 
peratures exceeding 24 K, selectively destroy HNC in favour of 



HCN. Interestingly, iGreaves & Nvmaiil ( [19961) find HCN/HNC 
abundance ratios of up to 6 for Galactic spiral arm clouds 
(GMCs). The highest ratios are observe d in regions of low - 
est N(H2), consistent with models by iSchiUce et alj (|1992|) . 
which show and increase of HCN/HNC with increasing C/CO. 
Bright HNC emission is commonly observed in extragalactic ob- 
iects (Aalto et al. 2002: . Wang et al . 2004; Meier & T urner 20051 ; 
iPerez-Beaupuits et^Tf 2007 ). In particular, overluminous HNC 
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7=3-2 is found in LIRGs (e.g., lAalto et al.ll2007h where gas 
tempe ratures, derived by the IR dust continuum ('Evans et al.' 
l2003h and mid-IR molecular absorption (Lahuisetal. 2007), 
are usually > 50 K and can reach values as high as a few 
100 K. Here, ion-molecule ch emistry in PDRs may b e respon- 
sible for the observed ratios. iMeiierink et al.l (l2007h find that 
the HNC/HCN 7=1-0 line ratio is enhanced in PDRs and can 
reach a maximum value of one for H2 column densities ex- 
ceeding 10^^ cm"^. Even higher HNC/ HCN fine ratios can re- 
sult from XDR emission. Models by iMeijerink et al.l (l2007h 
show that for the 7=4-3 transition, this ratio can be as high 
as 1.6 in dense (n(H)>10'' cm""*) X-ray-dominated regions. The 
average HNC/HCN 7=1-0 fine ratio for AGN and LIRGs in 
our sample is ^0.5 (s ee plot c in Fig . [Hi, w hich is consis- 
tent with the results by IMeiierink et alj (l2007h for XDRs with 
densities <10^ cm"^ and high X-ray fluxes. Fx >10 erg cm"^ 
s For the same sources, low-density XDR chemistry could 
also explain the observed HCO+/HCN ratios. The HNC/HCN 
ratio of order unity observed in UGC 5101 and NGC 2273 
(see Table |2]) is consistent with the emission coming from ei- 
ther a FDR or a low-density XDR (nn <10^ cm"^^). Models 
of PDRs with different density and radiation field strength can 
explain HNC/HCN 7=1-0 line ratios ra nging from 0.5 to 1 
(IMeiierink & Spaaiisll2005l; iLoenen et al. 112008). In our sample, 
the observed HNC/HCN ratio reaches values as low as 0.2 (plots 
b, c in Fig. |4]i, which cannot be explained by existing PDR or 
XDR models. However, line ratios HNC/HCN<0.5 are observed 
in Galactic PDRs. Abundance ratios of about 0.2 are found, e.g. 
, by Fuente et al. ( 1993) in the PDR region of the reflection neb- 
ula NGC7023. Low values o f HNC/HCN are also found in other 
galaxies (see , e.g., iBaan et al.. (2008,)) as discussed in Section 
1431 



5.4. Luminous HC3N: a nascent starburst tracer? 

Bright HC3N emission is often observed in Galactic ho t 
cores (iRodriguez-Franco et al.1 119981: Ide Vicente et al.1 12OOOI) . 
The molecule is easily destroyed by UV radiation and reactions 
with C^ and He^ (e.g.. Turner et al. 1998) and can only survive 
in regions shielded by large gas and dust col umns. In a high- 
resolu tion study of the chemistry of IC 342, iMeier & Turned 
(l200l find that HC3N emission follows the 3 mm contin- 
uum dust emission and anti-correlates with regions of in- 
tense UV radiation, su ch as PDRs. Observations in the Galaxy 
(IWvrowski et al.llT999h and in the LIRG NGC 4418 dAalto et all 
I2OO7L ICostagliola & Aaltd 12010) show that HC3N is strongly 
connected with the IR field by means of its vibrational bend- 
ing modes. In our sample, only five galaxies have detectable 
HC3N emission: NGC 4418, IRAS 17208, IC 860, NGC 1068, 
and NGC 7771. For completeness, we also include in our analy- 
sis the de tection of HC3N 7=10-9 in the ULIRG Arp 220 by 
lAalto et al. (2002). Since these objects span a wide range of 
LSR velocities, in most cases the detectability of HC3N does 
not depend on the sensitivity of our measurements, but on intrin- 
sic properties of the observed galaxies. This is further discussed 
in Appendix lAl where we compare HCN and HC3N detection 
thresholds. 

All galaxies with detected HC3N have HCO+/HCN <1 and 
higher-than-average HNC/HCN ra tios. This is consistent with 
the hot core models described in iBavet et al] (|2008|) . as dis- 
cussed in Sect. 15.21 It is thus possible that the bulk of the molec- 
ular emission in these objects comes from regions of young 
star-format i on. T his interpretation was already proposed by 
lAalto et al.l (l2007h to explain the HC3N emission of NGC 4418 



in terms of a nascent starburst. 

The HC3N-bright galaxies reside in the HCO+/HCN interval 
that in the diagnostic diagrams discuss ed by Kohno e t al. ( 20()l|) 
is associated with S eyfert galaxies. iKohno et al.l (l200lh and 
llmanishi et alj (l200 7i) interpret the observed ratios as the sig- 
nature of XDR chemistry driven by the hard radiation from 
the AGN. Our data seem to challenge this interpretation, since 
HC3N cannot survive in significant abundance in highly irradi- 
ated gas. Bright HC3N emission can emerge from an AGN only 
if the central source is shielded by a considerable amount of gas 
and dust. In these objects, the bulk of the molecular emission 
would be coming from an hot core-like region instead of be- 
ing dominated by XDR chemistry. Plots e and / in Fig. |4] also 
show that the galaxies with the highest HC3N/HCN ratios have 
high ratios of IRAS 60 to 100 fim fluxes. This may indicate that 
HC3N is more abundant in warm environments. 
An enhanced HC3N emission at high gas temperatures could 
also be owing to excitation effects. The upper state energy of 
the observed HC3N transitions is indeed of the order of 20-30 K, 
higher by a factor of ^5 than typical energies for other high den- 
sity tracers, meaning that the molecule is more efficiently excited 
at higher temperatures. The bright HC3N emission is thus likely 
to be a combination of enhanced abundance and favourable ex- 
citation conditions in warm, dense gas. 

Radiative pumping of the molecule's rotational levels by 
means of IR vibrationa l modes was also suggested by 
ICostagliola & Aaltol (|2010|) . 

5.5. COA^CO vs dust temperature 

Plot g in Fig. |5] shows how the ratio ^Ri_() =CO/'^CO 7=1-0 
varies with dust temperature. The dust temperature was calcu- 
lated from IRAS fluxes at 100 and 60 //m, via the formula 



Tdust = -(1 +z) 



82 



In (0.3/60;jm//l00//m) 



■0.5 



(2) 



see, e.g., Solomon et al.' (Il997h . The IRAS fluxes were obtained 
from the NASA/IPAC Extragalactic Database. 
A clear trend is observed, with !Ri_o increasing as Tdust 

increases. This is a well established resu lt, see, e.g., 

lYoung & SandersI (Il986h and lAalto et al.1 d 19951) for a detailed 
discussion. The trend is generally attributed to opacity effects 
arising in a collection of non-uniform clouds, averaged in the 
beam. In a molecular cloud, the emission of CO 1-0 arises 
mostly from the outer layers, up to optical depths t ^1, while 
the '^CO emission stays in general optically thin deeper into the 
volume of the cloud. The line ratio ''Xi_o is thus mostly affected 
by the relative opacity of the two species. 
The opacity of the CO 1-0 line strongly dep ends on the excita- 
tion temperature (see, e.g., ' Aalto et al.l 1995 1 Eq. 1). With rising 
Tex, the 7=1 level is in fact de-populated, causing a decrease in 
optical depth of the CO 1-0 transition. Where tco -1, we can 
assume rnco <<1 ™d the intensity of the '^CO 7=1-0 line is an 
even steeper function of Tex th an the '^CO intensity (see, e.g., 
iMauersber ger & Henkell ll'993h . As a consequence, ^Ri-o must 
increase with tempera ture. 

In our Galaxv iYoung & SandersI (11986! ) find that ^Ri-o has an av- 
erage value of roughly 8, whil e, in ext ragalactic surveys, values 
!Ri_o >20 are observed (e.g., lAalto et al. 1995, this work), the 
highest ratios being found in interacting systems and mergers. 
In our sample, the highest !Ri_o - 28 is found for the galaxies 
NGC 6240 and NGC 1614. In these objects, high gas tempera- 
tures and a diffuse ISM may result in an extremely low optical 
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depth of the CO 1-0 line, which would explain the observed ra- 
tios. 

Additional processes have been proposed in order to explain 
these extreme ratios. In the central regions of merging and in- 
teracting galaxies, large amounts of gas and dust are trapped 
in the deep potential well, giving rise to high gas pressures. 
Large turbulent motions, powered by differential rotation or stel- 
lar winds from an enhanced star-formation are required to sus- 
tain the clouds against collapse. This would cause a broadening 
of the CO emission line and a decrease in optical depth, with a 
consequent increase in ^Ri_o (.Aalto et al. 1995). 
Chemistry could also provide an explanation of the observed 
high ratios. In star-forming regions the surfaces of molecular 
clouds are indeed dominated by UV-photons from young stars. 
In these regions, '^CO is dissociated more easily than CO, which 
is less affec ted by the UV radiation thanks to self-shielding 
(lAalto etal.i ri995). This process could be relevant in starburst 
galaxies, as NGC 6240 and NGC 1614 in our sample. 
The age of the starburst could also have an influence on ^Ri_o. 
Nuclear processing in stars favours '^CO over CO, and the ev- 
idence for an e nhanced ''Xi-o at low me tallicities has been re- 
ported by, e.g., lYoung & SandersI (Il986l) . A large ''Xi_o would 
be thus compatible with a young star-formation age. 

5.6. CO/CNvsCO/C^^O 

Plot h in Fig. |5] shows the ratio CO/CN vs CO/C'**0. The total 
CN intensity was obtained by adding the integrated intensities 
of the two lines of the spin doublet 7=3/2-1/2 and 7=1/2-1/2. 
The ratio CO/CN varies in our sample over roughly one order 
of magnitude and shows a general trend of increasing CO/CN 
with increasing CO/C'**0. This implies that CN and C"^0 emis- 
sion are positively correlated. There are several potential expla- 
nations of this: 

- We can assume that both C'^'O and CN are optically thin, for 
CN we know this from the CN 1-0 spin group ratio, which 
is at the optically thin level for all detected galaxies (apart 
from, possibly, UGC 5101). In this case, the intensity of CN 
and C'^O would scale with the total molecular column den- 
sity. This linear scaling with total column would not hap- 
pen for the CO 1-0 intensity, which can be considered to be 
at least moderately optically thick in all the galaxies in our 



sample. The observed correlation would then be a result of 
optical depth effects. 
- It is also possible to speculate on other explanations. It is not 
immediately self-evident that the C'^O and CN column den- 
sities should follow each other. In PDRs, the CN molecule is 
usually abundant up to a visual extinction Av-2. At higher 
depths in the clouds the radical CN may react with other 
molecules to form - for example - HCN. Thus, CN can, to 
some degree, trace the amount of FDR surfaces. 
In PDRs, it is likely that C'**0 could become photo- 
dissociated. However, if the CN is tracing PDRs, it is also 
tracing the impact of the star-formation on host clouds. 
These are the regions where the ISM enrichment caused by 
the starburst hits first., i.e., where C'**0 would become en- 
hanced early on. Hence, in this scenario, a C'^O-CN corre- 
lation would be expected as a result of ISM enrichment from 
the starburst. 

In this context, it is interesting to note the distinct deviations 
from this correlation. We have the three outliers with faint C'^O 
compared with their CN luminosities: NGC 4194, NGC 1614 
and NGC 7469. A possible way to obtain this comparable over- 
luminosity in CN is to have CN emerging from an AGN compo- 
nent, with no associated enrichment from a starburst component. 
The C'^O enrichment should come quickly in a starburst. If '^O 
is a product of a-capture on '"^N during He burning, it will come 
from mass loss in He-burning systems, such as Wolf-Rayet stars 
and SNe type II. This means that it should come up with prompt 
nucleosynthesis - e.g., follow O abundance. 
Because the molecular emission from NGC 7469 is dominated 
by the starburst ring, this comparatively C'^O weakness is some- 
what puzzling and should be further investigated. The observed 
high CO/C"^0 ratio could also be caused by a low optical depth 
of the CO 1-0 line. 

The other two outliers are different since they are both merg- 
ers where the CO emission is dominated by a prominent cross- 
in g dust lane and no t the starburst region (Olsson et al. 201Q|; 
lAalto & Hiittemeistej l2000h . It is therefore possible that the 
molecular gas in the IRAM beam is dominated by in-falling gas 
in the d ust lane, which has n ot yet been enriched by the starburst. 
Indeed. 'Marti n et al.1 (1201 Ol) find high '^C/C ratios from CO and 
CCH observations in the starbursts M 82 and NGC 253, suggest- 
ing that the bulk of the molecular mass in these galaxies mostly 
consists of unprocessed material. 
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An outrigger in a different direction is IC 860. This galaxy has 
the lowest CO/C'^0 line ratio in the sample - but the highest 
CO/CN 1-0 line ratio. In this source, it is possible that only a 
small fraction of the cloud surface area is in the form of PDRs - 
and/or that the CN is being used up in forming HC3N together 
with C2H2, which is boiling off the grains. Its luminous PAH 
emission (see Fig. |4j/) does suggest however that PDRs should 
be quite abundant. It is possible that the PDRs in this object do 
not coexist with regions dense enough to excite CN emission. 

5. 7. Other isotopic line ratios 

For a small subset of galaxies, it was possible to obtain isotopic 
line ratios for HCN and NCO^, which are reported in Table |5] 
If we assume a Galactic C/'"'C abundance ratio of ^50 (see, 
e.g., IWilson & Rood 1994), or as high as 100 (Martin et al. 
l2010l) . it is evident that in the galaxies where '^C variants were 
detected, HCN and HCO^ 7=1-0 transitions are at least mod- 
erately optically thick, with HC*/H'^C* as low as 4. These 
data show how the assumption of optically thin lines when in- 
terpreting extragalactic molecular intensity ratios is not always 
correct and opacity effects have to be included in the analysis, 
even for species much less abundant than CO. The broad band 
of EMIR allows us to determine isotopic intensity ratios with 
unprecedented accuracy, which will be of great importance in 
constraining models of galactic emission. However, deeper inte- 
grations are still needed to set reliable limits in more optically 
thin environments. We also advise caution when interpreting 
HCOVH'^CO+ intensity ratios, since the H'^CO+ 7=1-0 Une 
is blended with SiO emission and the two lines are not easy to 
separate at low SNRs. 



Galaxy 


HCN/H'^CN 


HCO+/H'^( 


NGC 1614 


>2.8 


4(2) 


Mrk231 


6(2) 


7(6) 


NGC 4418 


8(3) 


>4 


NGC 1068 


17(7) 


24(20) 



Table 5. Isotopic line intensity ratios for HCN and NCO^. 
Uncertainties (1-cr) are shown in parenthesis. 



measurement errors. The 7=10-9 V7=l/ transition at 91 GHz 
has no detected e counterpart and will therefore be treated as an 
upper limit in our excitation analysis. 

A first analysis of the excitation of the molecule was performed 
via the population diagram method. This results in an excitation 
temperature of 15 K for the v=0 transitions and 42 K for the 
V7=l lines. 

A rotational te mperature of 15 K is significantly lower than the 
30 K derived bv lCostaghola & Aa lto (2010) for the other HC3N- 
bright LIRG, NGC 4418. This may suggest that bright HC3N 
emission could originate from different environments. However, 
our estimate is based only on two transitions and has to be fur- 
ther investigated. 

The vibrational excitation of HC3N is strongly dependent on the 
IR field, since the molecule's bending modes can only be ex- 
cited by radiation. A higher excitation of the V7=l lines, com- 
pared with the v=0 levels, is therefore expected in a source with 
a strong IR emission such as IC 860. The derived 42 K is a lower 
limit to the excitation temperature, because we only have an up- 
per limit for the 7=10-9 v-j-lf transition at 91 GHz. 
In highly obscured galaxies, the detected IR emission is in most 
of the cases optically thick. Because millimeter radiation can 
penetrate a larger column of gas and dust, vibrationally excited 
HC3N is a valuable probe of the IR field beyond its optically 
thick surface. 

The first multi-transition analy sis of extragalactic vibration ally 
excited HC3N was reported by ICostaghola & Aaltol ( 1201 Oh for 
the IR galaxy NGC 4418. Here the authors find a HC3N vibra- 
tional temperature of 500 K, mu ch higher than the temperature 
of the optically thick dust (80 K, Evans et alj (l2003h ). which may 
be the evidence of an obscured compact object in the centre of 
the galaxy. 

It would be interesting to perform this kind of analysis on IC 860, 
which presents unresolved emission in both IR and radio and 
would then be a good candidate for a deeply embedded compact 
source. 



Transition v jT^dv Av E^p 

[GHz] [Kkms-'] [ km s"' ] [K] 



5.8. HCsNin IC 860 

In total we detected four emission lines of HC3N in IC 860, re- 
ported in Tabled In the 112 GHz band, we detected two tran- 
sitions of the vibrationally excited state V7=l (e,/), together 
with the 7=12-11, v=0 Une. In the 88 GHz band we detected 
the 7=10-9, v=0 line, and we have a tentative detection of the 
V7=l/, vibrationally excited transition. 

The interaction between the bending angular momentum of the 
vibrationally excited states and the rotational angular momen- 
tum of the molecule leads to a 1-splitting of the levels. Each vi- 
brational state is thus split into two levels, labelled e or / depend- 
ing on the wavefunctions parity properties. These parity labels 
are also shown in Table |6] A useful r eference for the labe lling of 
doubled levels in linear molecules is lBrown et al. l (fT97l . 
Both e and / states have similar excitation properties and should 
result in an emission doublet with comparable line intensities. In 
our spectra, this is true for the 7= 12- 11 V7=le and V7=l/ lines 
at 109 GHz, which show integrated line intensities within the 



J= 


10-9 




90.978 


0.43(0.09) 


200 


24 


7= 


10-9 V7= 


1/ 


91.333 


<0.12 


170 


342 


7= 


12-11 




109.173 


0.75(0.16) 


300 


34 


7= 


12-11 V7 


= le 


109.441 


0.33(0.08) 


170 


355 


7= 


12-11 V7 


= 1/ 


109.598 


0.21(0.08) 


170 


355 



Table 6. Detected lines of HC3N in IC 860. The integrated in- 
tensity is not corrected for main beam efficiency (see Table[T]i. 



5.9. Models of galactic molecular emission 

The line ratios observed by iBaan et al ] (120081) were co mpare d 
with a grid of XDR and PDR models by iLoenen etal] (|2008|) . 
These authors find that the low HNC/HCN<0.5 line ratios can- 
not be explained by XDR or PDR chemistry, and they claim that 
a possible explanation could come from mechanical heating. At 
temperatures > 30 K, neutral-neutral reactions efficiently destroy 
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HNC in favour of HCN (iHirota et al.ll998l) . lLoenen et alJ (l2Q08h 

suggest that shocks from supernova explosions may heat the gas 
to temperatures >100 K, resulting in a depletion of HNC. The 
observed HNC/HCN line ratios can be explained by mechanical 
heating for a gas density of 10**^ and a star-formation rate (SFR) 
of 20 Mq yr^'. Higher densities would require higher SFR, 
which would result in the destruction o f the mo lecular material 
surrounding the star-formation region dLoene n et al. 2008). In 
this scenario, the observed low HNC/HCN line ratio would orig- 
inate from low-density («// <10^ cm"-'), warm (T^lOO K) gas. 
The corresponding high HCO^/HCN observed in plot Fig. HJ? 
may be a consequence of the low gas density. As discussed in 
Sect.ES HCO+ 1-0 has a lower critical density (^10** cm"^^) 
compared with HCN 1-0 (^10^ cm~^) and is more efficiently 
excited in diffuse gas. An increased HCO^ abundance may also 
be caused by dissi pation of turbu l ence driven by the supernovae 
winds. Models by iJoulain et alj (1199 8) show that in low den- 
sity gas, bursts of viscous dissipation can heat the gas up to 
the high temperatures required by endothermic reactions, such 
as C+-HH2 ^CH+-hH {hElk ^4600 K), more efficient than the 
standard cosmic ray-driven chemistry. This leads to an enhance- 
ment of CH^, OH, and HCO^, and may help to create the high 
HCO^/HCN ratios we observe. 

This picture is comp licated by models of HCN and HNC chem- 
istry as discussed in ISchiUce et al.l (Il992h . For gas temperatures 
- 100 K, an efficient destruction of HNC by neutral-neutral reac- 
tions requires densities higher than 10^ cm while at lower den- 
sities the HCN/HN C abundance ratio is about unity (Fig. 13 in 
ISchilke et al.ll 19921) . However, at higher gas temperatures (>200 
K), HNC destruction may become effective even for densities as 
low as ^10"* cm"^ and thus provide an explanation for the low 
HNC/HCN. 

5.10. HCCr, HNC, and PAHs 

In Fig. lU^, we compare the observed HNC/HCO^ 7-1-0 line 
ratios with the PAH equivalent widths (EW) from (ISpoon et al.l 
12007 ). The distribution of HNC/HCO+ can be described as a step 
function: it is roughly constant for low PAH EWs and has a sud- 
den drop for LogioEW>-0.5. The galaxies with high PAH EWs 
are mostly starbursts or LIRGS. 

PAHs are excited by UV radiation and destroyed by hard radia- 
tion (e.g.. X-rays). Because if this, strong PAH emission is usu- 
ally interpreted as a clear signature of star-formation, instead of 
AGN activity. A decrement of HNC/HCO+ for high PAH EWs 
may then imply either that HCO^ is enhanced in a PDR, or that 
HNC is depleted. Indeed, if we compare the HCO+ and HNC 
emission to HCN 1-0, for increasing PAH EWs we have an in- 
crease in the HCO^/HCN ratio and a decrease in HNC/HCN. 
This can be seen in the graphs in Fig. IB. II 
If directly connected to the molecular abundances, the observed 
ratios are difficult to explain with P DR and XDR chemistry. 
Models bv Meiier ink & SpaansI 12005) show that HCO+/HCN is 
lower in PDRs than in XDRs, and HNC/HCN very close to unity 
in a wide range of densities and gas temperatures. Howev er, the 
observed trend agrees with results by Imanishi e t alJ (12007 '). who 
find that AGNs have a lower HCO+/HCN ratio than starburst 
galaxies. 

A depletion of HNC is also expected in shocks. Here, because 
of the high temperatures and densities, neutral -neutral reac- 
tions efficiently destroy HNC in favour of HCN ( SchiUce et alJ 
Il992h . Mechanical heating from shocks generated by super- 
nova explosions may explain the observed trend. In this case, 
the HNC/HCN ratio should follow the dependence on temper- 



ature of neutral-neutral reactions. In Fig. we compare the 
HNC/HCN ratio with the ratio of the IRAS fluxes at 60 and 
100 yum. The graph clearly shows that there is no strong de- 
pendence of the HNC/HCN ratio with dust temperature. At 
densities >10^ cm""*, gas and dust thermalise and this should 
show in the diagram, since neutral-neutral reactions strongly 
depend on gas temperature. At densities below 10^ cm""*, the 
thermal coupling becomes poor, however, and it is possible 
to have Tg^s >Tdust- This scenario is complica ted by observa- 
tions of HCO^ in shocked gas in the Galaxy (|j0rgensen et al.l 
2004), which reveal that the molecule is destroyed by reactions 
with H2O in the shock front. Detectable HCO^ abundances are 
found i n the aftermath of the shock, where water freezes on 
grains. |j0rgensen et all ( |2004|) report an HCO^/HCN column 
density ratio of 0.6 in the outflow of the low-mass protostar 
NGC 1333 IRAS 2A. 

6. Conclusions 

We observed the molecular emission of 23 galaxies with the 
new EMIR receiver at IRAM 30m in the frequency bands 88 
GHz and 112 GHz. The observed line ratios were compared 
with existing models of emission from PDRs, XDRs, and 
hot cores. All non-compact starburst galaxies are found to be 
HCO^ -luminous, with Hne ratios HCO^/HCN > 1. Luminous 
infrared galaxies and Seyferts do not show any strong trend. 
However, the majority of active galactic nuclei do reside at low 
HCO^/HCN values compared with starbursts. 
When HNC/HCO^ increases, galaxies move towards higher 
HNC/HCN and HCO+/HCN ratios. Galaxies with HCO+/HCN 
< 1 present on average a 30 % increase in HNC/HCN, com- 
pared with those at higher HCO^/HCN ratios. If we exclude 
die Seyfert NGC 2273, the brightest HNC emission is found 
in luminous infrared galaxies. Starbursts are generally charac- 
terized by faint HNC and bright HCO^ emission (compared 
with HCN), while LIRGs mostly occupy the opposite end of the 
HNC-HCO+ con-elation, with high HNC and low HCO+ inten- 
sities. All HNC -bright galaxies are either LIRGs or Seyferts. 
AU HC3N detections have HCO+/HCN < 1 and HNC/HCN>0.4. 
Most of these are are LIRGs, the only exception being the 
Seyfert NGC 1068. 

Starburst galaxies tend to have high HCO+ /HCN and high PAH 
EWs, in agreement with observations by iKohno et~aL (1200 ll) 
and Imanishi et ak (2007,) . 

The interpretatio n of the low HCO^/H CN ratio as an effect of 
XDR chemistry ('iman ishi et alJl2007h is inconsistent with our 
observations, since we find bright HC3N emission in galaxies 
with the lowest HCO+/HCN. The HC3N molecule is easily 
destroyed by reactions with ions and hard UV and is thus 
extremely u nhkely to su r vive in X-ray dominated environments. 
Models by 'B avet et alj (12008*) for hot cores provide a better 
explanation for ffie observed ratios, suggesting that bright HC3N 
and faint HCO^ emission might be good tracers of embedded 
star-formation. 

Vibrationally excited HC3N was detected in IC 860, where it 
was possible to obtain an estimate of the rotational temperature 
of the v=0 transitions of 15 K, and an upper limit for the vj=l 
levels of 42 K. This galaxy has the highest HC3N/HCN ratio 
of the sample, representing an interesting target for further 
investigation of the molecule's properties. 
In agreement with previous studies, we find that the CO/'^CO 
intensity ratio is positively correlated with dust temperature, 
with typical values ^ 20, higher than the ones observed in the 
Galaxy. We suggest that these high values are caused by either 
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high gas temperatures or infall of unprocessed material in the 
centres of starbursts. 

The emission of CN and C'^O appear to be positively correlated. 
This may simply be because both CN and C^^O are optically 
thin and scale with the total gas column, or may be caused 
by ISM enrichment by the starburst activity. CN emission is, 
indeed, particularly strong in PDR surfaces, which are supposed 
to be tracing the impact of star-formation on the host clouds. 
In agreement with previous extragalactic studies (e.g., 
lBaan etal.1 l2008h . we find HNC/HCN intensity ratios lower 
than 0.4, which cannot be explained by current models of PDR 
and XDR che mistry. A possible explanation was suggested by 
iLoenen et al.l (2008) to be mechanical heating by supernova 
explosions. This model would require low gas densities (<10^ 
cm~^) or extremely high star-formation rates (>20 M©), but 
would explain both the low HNC/HCN and high HCO+/HCN 
ratios in starbursts. Mechanical heating by SN explosions could 
also explain the observed drop of the HNC/HCO^ ratio for 
R\H EWs > 0.3 fim, if we assume that PAHs efficiently trace 
starformation. 

However, a temperature dependence of the HNC/HCN ratio, 
which would be expected if the HNC abundance was depleted 
by neutral-neutral reaction in the shocked gas, does not emerge 
from our data. Furthermore, it is not yet clear how SN shocks af- 
fect HCO^ abundance, with diff'erent models giving contrasting 
results. 

6. 1 . Future directions 

Evidently, explaining the molecular emission observed in galax- 
ies by means of the existing chemical models is a challenging 
task. Our understanding of the gas physical conditions is still af- 
fected by many uncertainties. Mechanical heating, invoked to 
explain the low HNC/HCN ratios, requires further investiga- 
tion with shock models and observations of shock tracers (e.g., 
IR transitions of H2 or SiO). For single dish observations, the 
analysis of the observed emission is complicated by the lack 
of spatial information. Radiation from different molecules and 
transitions may in fact be coming from different regions in the 
galaxy, and thus it may be impossible to explain observations 
with only one set of physical conditions. Spatial distribution and 
density gradients in the molecular gas may have a strong impact 
on the observed line ratios. This may be particularly important 
for the HCO^/HCN ratio, which appears to be strongly density- 
dependent, because of both chemical and excitation effects. The 
angular resolution needed to resolve the spatial distribution of 
molecular gas in distant galaxies will become available in the 
foreseeable future with the ALMA interferometer Observations 
of mm lines of formaldehyde (H2CO) may also help to con- 
strain density and temperature in these objects. Sub-mm obser- 
vations of high-J transitions of CO, which have been indicated 
as excellent tracers of PDRs by Meiierink et al. (2007), are now 
available with the Herschel telescope. The increased sensitiv- 
ity of new generation mm- and sub-mm observatories will al- 
low us to detect less abundant molecular species (e.g., HNCO, 
CH3OH) that may be better tracers of interstellar chemistry (e.g., 
[Meier & Turner 2005; Martin et al. 2009). The evolution of the 
starburst phase can be studied with radio continuum observa- 
tions, which will be discussed in Paper II. We also plan to further 
explore molecular excitation in our sample by means of multi- 
transition studies. We successfully applied for telescope time at 
IRAM 30m to observe 1 mm lines of the species discussed in 
the present work. Collected data will be discussed in Paper III, 
where we will present a NLTE analysis of the emission of the 



main molecular species, together with a comparison with a grid 
of PDR and XDR models. 
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Fig. A.l. Ratio of HC3N 7=10-9 over HCN 7=1-0 integrated 
intensities against the HCN 7=1-0 signal-to-noise ratio. The line 
shows a linear fit through the HC3N detections. See Fig.|4]for an 
explanation of the symbols. 
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Appendix A: Detection of I-IC3N 

The ratio of HC3N 7=10-9 over HCN 7=1-0 integrated intensi- 
ties is plotted in Fig. lA.ll against the HCN 7=1-0 signal to noise 
ratio (SNR). The graph clearly shows that the non-detections are 
not caused by the sensitivity limit of our observations. In a wide 
range of SNR values, detections and non-detections coexist. 
The line in Fig. lA.ll shows a first order fit through the HC3N de- 
tections, excluding the low-SNR detection in NGC 7771 (num- 
ber 12 on the graph). Our detection limits (arrows in the graph) 
are, on average, 3 times fainter than the expected intensity de- 
rived by the fit. This means that in most of the cases - this may 
not be true for NGC 2273, number 19 in the graph - we do not 
detect HC3N not because of a lack of sensitivity, but because the 
emission is intrinsically faint. In NGC 1068 the molecular emis- 
sion is about one order of magnitude brighter than in the other 
galaxies of the sample, therefore the HC3N detection (with the 
lowest HC3N/HCN ratio in the sample) in this galaxy may be 
only caused by its intrinsic high luminosity. 

Appendix B: HCO , HNC and PAHs 

The strong dependence of the HCO+/HNC ratio with PAH 
equivalent width, shown in Fig. |4j/, can be better understood by 
plotting the behaviour of the single tracers with increasing PAH 
EW. 

In Fig. 10 we plot the ratios HCO+/HCN and HNC/HCN 
against the PAH strength. Evidently, at high PAH EWs we have 
a depletion of HNC (respect to HCN) and an enhancement of 
HCO^. This two trends both contribute to the fast decrease of 
the HCO^/HNC ratio seen in Fig.|4j/. For a discussion, see text 
in Section lSTTOl 



Log.i(j PAH equivalent width (6.2 urn) [nm] 

Fig.B.l. Line ratios of HCO+ and HNC over HCN against the 
PAH equivalent width. Notice the decrease of HNC and increase 
of HCO+ at high PAH EWs, to be compared with Fig.|4j/. 

Appendix C: Comparison with the literature 

In Ta blelC.l | we compa re our observations with literature data 
from iBaan et aP (l2008h . Antenna temperatures, obtained with 
EMIR, were translated into flux densities, by applying the Jy/K 
conversion factor 

J[Jy] = niK] X 3.906 X Feff/77A, (C.l) 

with Feff forward efliciency, and t/a antenna efficiency of the 
30m telescope at 88 GHz. 

Appendix D: ElVIIR Spectra 

In Fig. ID.ll|D.9| we show the spectra observed in all the galax- 
ies in the sample. The intensities are in T^, not corrected for 
beam efficiency (see Table[TJ. On the x-axis we plot the sky fre- 
quency in GHz, not corrected for redshift. The position of the 
main emission lines is also shown. 
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Table C.l. Comparison with data from the hterature, as reported by Baan et aP (l2008h (BOS), and this work (EMIR). Exact refer- 
ences for the observations can be found in Table Bl in lBaan et al. I (12001 7 Differences (in %) between BOS and EMIR line intensity 
ratios are shown in column A%. 
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Fig.D.l. Observed spectra in the 88 GHz band. The intensity scale is in T^, not corrected for main beam efficiency. The main 
molecular transitions are labelled regardless of line detection. The CiH 3/2 and 1/2 labels mark the limits of the CiH multiplet at 87 
GHz. Transitions of vibrationally excited HC3N are labelled as v^le and V7I/. The frequency scale is the observed frequency, not 
corrected for redshift. 
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Sky Frequency [GHz] 
Fig.D.2. Observed spectra in the 88 GHz band (continued). See caption of Fig. ID. II 
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Fig.D.3. Observed spectra in the 88 GHz band (continued). See caption of Fig. ID. II 
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Fig.D.4. Observed spectra in the 88 GHz band (continued). See caption of Fig. ID. II 
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Fig.D.5. Observed spectra in the 88 GHz band (continued). See caption of Fig. lDTI 
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.D.6. Non-detections in the 88 GHz band. See caption of Fig. ID. II 
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Fig.D.7. Observed spectra in the 112 GHz band. The intensity scale is in T^, not corrected for main beam efficiency. The region 
marked with (1/10) around the CO 1-0 line was scaled down by a factor 10. The main molecular transitions are labelled regardless 
of line detection. Transitions of vibrationaUy excited HC3N are labelled as v^le and V7I/. The frequency scale is the observed 
frequency, not corrected for redshift. 
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Fig. D.8. Observed spectra in the 112 GHz band (continued). See caption of Fig. ID. 71 
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Fig. D.9. Observed spectra in the 112 GHz band (continued). See caption of Fig. ID. 71 
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Table E.l. Summary of the line properties of the observed galaxies. Integrated intensities and line widths were derived by means 
of Gaussian fitting. Uncertainties are shown in parenthesis. (*): For C2H, the integrated intensity is obtained by fitting a single 
Gaussian to the multiplet at 87 GHz. This results in large line widths, which are not to be attributed to dynamics. 
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Table E.2. Summary of the line properties of the observed galaxies (continued). See caption of Table [E. II 
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<0.51 






Si02-1 


<0.18 


- 


HC3N12-11 


<0.51 


- 




H13CN1-0 


<0.18 




HC3N12-llv7e 


<0.51 






C2H1-0 


<0.18 




C18O1-0 


0.51(0.10) 


238(50) 




HCNl-0 


<0.18 




13CO1-0 


1.40(0.10) 


196(13) 




HCO+1-0 


<0.18 




CNl-0J=l/2 


0.51(0.14) 


238(62) 




HNCl-0 


<0.18 




CNl-0J=3/2 


1.10(0.20) 


317(67) 




HC3N10-9 


<0.18 




CO 1-0 


20.25(0.20) 


200(16) 


UGC2866 


H13CO+1-0 


<0.21 




CH30H 


<0.60 






Si02-1 


<0.21 


- 


HC3N12-11 


<0.60 


- 




H13CN1-0 


<0.21 




HC3N12-llv7e 


<0.60 






C2H1-0 


1.30(0.15) 


397(50)* 


C18O1-0 


0.60(0.14) 


283(56) 




HCNl-0 


1.47(0.12) 


254(25) 


13CO1-0 


3.46(0.12) 


283(10) 




HCO+1-0 


2.15(0.12) 


273(15) 


CNl-0J=l/2 


0.69(0.12) 


190(30) 




HNCl-0 


0.75(0.13) 


336(55) 


CNl-0J=3/2 


1.98(0.16) 


264(22) 




HC3N10-9 


<0.21 




CO 1-0 


71.54(0.28) 


266(1) 


UGC5101 


H13CO+1-0 


<0.31 




CH30H 


<1.85 






Si02-1 


<0.31 


- 


HC3N12-11 


<1.85 


- 




H13CN1-0 


0.37(0.17) 


456(130) 


HC3N12-llv7e 


<1.85 






C2H1-0 


0.98(0.27) 


1043(130)* 


C18O1-0 


<1.85 






HCNl-0 


1.10(0.15) 


455(130) 


13CO1-0 


<1.85 






HCO+1-0 


0.40(0.16) 


500(130) 


CNl-0J=l/2 


1.85(0.19) 


578(112) 




HNCl-0 


0.90(0.18) 


476(130) 


CNl-0J=3/2 


1.85(0.19) 


578(112) 




HC3N10-9 


<0.31 




CO 1-0 


16.20(0.30) 


536(112) 


NGC2273 


H13CO+1-0 


<0.33 












Si02-1 


<0.33 


- 










H13CN1-0 


<0.33 












C2H1-0 


<0.33 












HCNl-0 


0.55(0.10) 


154(66) 










HCO+1-0 


0.33(0.10) 


154(66) 










HNCl-0 


0.33(0.10) 


154(66) 










HC3N10-9 


<0.33 











Table E.3. Summary of the line properties of the observed galaxies (continued). See caption of Table IE. 1 1 
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Galaxy 


Transition 


Kkms"' 


km s"' 


Transition 


jT*dv Ay 
Kkms ' s ' 




88 GHz Band 






112 GHz Band 




NGC1377 




<0.3 


70 






IRAS 15250 




<2.7 


170 







Table E.4. Upper limits (3-cr) for non-detections, estimated from noise rms. The assumed line width was derived from CO 1-0 
observations. 



